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Quantum electrodynamics describes the interactions of elec- 
trons and photons. Electric charge (the gauge coupling constant) 
is energy dependent, and there is a previous claim that charge is 
affected by gravity (described by general relativity) with the impli- 
cation that the charge is reduced at high energies. But that claim 
has been very controversial with the situation inconclusive. Here I 
report an analysis (free from earlier controversies) demonstrating 
that that quantum gravity corrections to quantum electrodynamics 
have a quadratic energy dependence that result in the reduction of 
the electric charge at high energies, a result known as asymptotic 
freedom. 

The standard model of particle physics is based in part on quantised 
Yang-Mills (or non-Abelian) gauge fields. The quantisation of such fields 
leads to the important prediction that at high-energy the strength of the 
gauge coupling constant governing the interaction of the Yang-Mills fields 
weakens, a phenomenon known as asymptotic freedom-^. This contrasts 
with a theory like quantum electrodynamics (that describes electrons and 
photons) where the gauge coupling constant, the electric charge in this case, 
gets stronger as the energy scale increases. A practical consequence of this 
is that for Yang-Mills theories that are asymptotically free a perturbative 
approach based on a weak coupling constant becomes more reliable at high 
energy, whereas for quantum electrodynamics perturbation theory breaks 
down. 

The key equations that govern the behaviour of the coupling constants in 
quantum field theory are the renormalisation group Callan-Symanzik equa- 
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tions 3,4 . If we let g denote a generic coupling constant, then the value of g 
at energy scale E, the running coupling constant g(E), is determined by 

E*§- = «*,,), (1) 

where j3(E,g) is the renormalisation group /3-function. Asymptotic freedom 
is signalled by g(E) — >■ as E — > oo, requiring (3 < in this limit. 

In the standard model of particle physics gravity is usually ignored as 
it plays an inessential role in most calculations of interest. Additionally, if 
we view Einstein's theory of gravity as a fundamental theory it exhibits the 
undesirable property of non-renormalisability-~-, and hence lacks the power 
of predictability. However, it is possible instead to view Einstein gravity 
as an effective theory that is only valid below some high energy scale 1 ^ 1 -. 
The cut-off scale is usually associated with the Planck scale Ep ~ 10 19 GeV, 
the natural quantum scale for gravity. Above this energy scale some theory 
of gravity other than Einstein's theory applies (string theory for example); 
below this energy scale we can deal with quantised Einstein gravity and 
obtain reliable predictions^— . Adopting this effective field theory viewpoint 
it is perfectly reasonable to include Einstein gravity with the standard model 
of particle physics and to examine its possible consequences using quantum 
field theory methods. 

With this in mind, a calculation was performed— in Einstein- Yang-Mills 
theory that looked at the effects of quantum gravity on the running Yang- 
Mills gauge coupling constant. This calculation found that quantum gravity 
leads to a correction to the renormalisation group /3-function (not present in 
the absence of gravity) that tends to result in asymptotic freedom; this holds 
even for theories (like quantum electrodynamics) that are not asymptotically 



free when gravity is neglected. Possible phenomeno logical consequences of 
this result have been considered—. 

The potential importance of the original calculation— stimulated a num- 
ber of further investigations that cast doubt on its findings. It was shown— 
that a different choice of gauge condition led to the absence of any quan- 
tum gravity correction to the Yang-Mills /3-function. Because of the pos- 
sible gauge condition dependence, a gauge-condition independent version of 
the background-field method— 1 ^ along with dimensional regularisation— was 
used— and also found no quantum gravity contribution to the gauge coupling 
constant. The original result-^ came from quadratic divergences that auto- 
matically get regulated to zero using dimensional regularisation—. 

The situation was analysed with a traditional Feynman diagram approach 
using standard Feynman rules and it was shown— that if a momentum space 
cut-off was used the quadratic divergences could be made to cancel leaving a 
result that was consistent with dimensional regularisation. On the other hand 
a novel method of regularisation, called loop regularisation, claimed— that 
the quadratic divergences were present, contradicting earlier results^ 1 ^ 1 ^. 
Subsequent work has looked at the role of the cosmological constant— 1 ^, 
corrections in higher dimensions^, and scalar— and Yukawa fields^. Ad- 
ditionally, related calculations have been performed in the exact renormali- 
sation group approach— - —. 

The issue of quadratic divergences and the possible gauge condition de- 
pendence is the central theme the present paper. The generalised background- 
field method—*^ that is used here differs from the usual one by the addition 
of an extra contribution that is essential for maintaining gauge condition 
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invariance. The difference between the standard and generalised effective 
actions can be controlled by a single parameter v that when set equal to one 
gives the gauge condition independent result, and when set to zero gives the 
standard result that coincides with a traditional Feynman diagram analysis. 
Gauge invariance is maintained by using a proper 

time 30 1 3i cut _ ff_ The im- 
portance of this choice of regularisation is that it is capable of revealing any 
possible quadratic divergences unlike many other methods. It will be shown 
that the quadratic divergences that are obtained do not depend on the pa- 
rameter v, and that the traditional background-field method leads to a gauge 
condition dependent result for the charge renormalisation. The contribution 
of the extra terms needed to obtain the gauge condition independent result 
is greatly simplified— by choosing a particular gauge condition described be- 
low. Because the generalised effective action is independent of the choice 
made for the gauge condition there is no generality lost by this procedure. 
(A comprehensive review and justification of the method can be found in 
a recent monograph 33 .) Surprisingly, it will be shown that the quadratic 
divergences do not cancel, and that the basic prediction originally found 13 
is correct, although with a slight modification of the /3-function. Quantum 
electrodynamics is found to be asymptotically free if it is coupled to quantum 
gravity. 

Effective action 

The model studied is Einstein gravity with a cosmological constant coupled 
to quantum electrodynamics in four spacetime dimensions. (The presence 
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of the cosmological constant is not essential to the point that will be made, 
but is included for generality.) The basic fields are the spacetime metric 
g^u, the electromagnetic field A^, and the Dirac field ip. The principle aim 
is to calculate the quantum gravity contribution to the renormalisation of 
the charge. To do this it is sufficient to adopt the background-field method, 
choose the background metric to be flat, the background Dirac field to van- 
ish, the background electromagnetic field to correspond to a constant field 
strength F^, and to concentrate on that part of the one- loop effective action 
that is divergent and involves the square of the background field strength. 

A standard calculation shows that the effective action to one-loop order 
is given by 

= - In det A* 3 - - lndetQ^ - In det (17^ + erfA^ - im). (2) 

The last term (with A^ the background gauge field) is the result of performing 
a functional integral over the Dirac field. The middle term is the contribution 
from the ghost fields required to remove the unphysical degrees of freedom 
of the gravity and electromagnetic fields. The first term is the result of 
integrating over the spacetime metric and electromagnetic fields; A l j is a 
second order differential operator that can be found from earlier work— & 
and will not be written down here due to its complexity. It is found by 
expanding the classical Einstein-Maxwell action about the background fields 
using g^ v = S^ + nh^ and A^ = A^ + a^. A Riemannian metric and standard 
conventions^ for the Dirac fields are chosen, k, 2 = 32ttG with G Newton's 
constant. 

The gauge-condition independent background-field method— 1 ^ is used as 
described in the final paragraph of the introduction. The gauge conditions 

6 



adopted are (with h = S^h^) 

Xx = d»h Xu -^d x h + ^u(A x d»a, + F, x a»), (3) 

X = d^. (4) 

To obtain the gauge condition independent result it is essential to take u — 1 
in our calculation. However u is kept general at this stage to indicate the 
gauge condition dependence of standard methods of calculation. The choice 
u = is often made and is called the de Donder gauge. 

The gauge conditions are incorporated in the usual way by altering the 
action with a gauge fixing term 

S gf = ^J d 4 x X X Xx + ^J^ X X 2 - (5) 
Here £ and ( are two dimensionless parameters. The choice £ = 1 = ( 
is often made because it simplifies the calculation enormously. However in 
doing this you lose any hope of addressing the gauge condition dependence 
of the calculation. £ and ( will be kept general here noting that the gauge 
condition independent effective action is the same as what is found by taking 
the limits £ — > and ( — > (along with u = v = 1). 



Heat kernel 

For operator A l j the heat kernel K t j(x,x'; r) is defined by 

- ^-ICjix^r) = AV^(z,x';t), ( 6 ) 

with boundary condition K l j(x, x'\ r = 0) = Sj5(x, x'). r is called the proper 
■k UTie 30 1 3i_ rp^g Q reen f unc tion G l j(x,x') for the operator A l j is 

A i k G k j(x,x') = Sj5(x,x'). (7) 
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It follows that the Green function and heat kernel are related by 

oo 



The importance of the heat kernel for quantum field theory arises from 
the existence^i^ of an asymptotic expansion as r — > 0: 

oo 

lCj{x, x; t) ~ (47rr)~ n/2 ^ T r E r l j{x) (9) 

r=0 

where n is the spacetime dimension (chosen as 4 here) and the heat kernel 
coefficients E r l j(x) depend only locally on the details of coefficients of the 
differential operator A l j. For many operators of interest in physics the results 
of the first few coefficients are known— however, the operators needed for 
the present calculation are more general than considered so far. (Some checks 
on our results follow from a different method—.) 

The divergent part of the effective action, as well as the Green function, 
can be related to the heat kernel coefficients. Formally (before regularisation) 

oo 

L A = ^]ndetA i j = ~J d n x /"^triT^^r). (10) 

o 

The one-loop effective action (T2J) is then given by 

= L A - 2L Q - 2L Dirac . (11) 

As with the Green function (j7j) the divergent part of (1101) comes from the 
t ~ limit of the proper time integral. To deal with this in a way that 
respects gauge invariance, general coordinate invariance, and the structure 
of any quadratic divergences, a proper time cut-off r c is used where the lower 
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limit of the proper time integration is replaced with r c . The divergences will 
show up 30 as r c — > 0. Because of this the divergent part divp(LA) of La 
follows from using the asymptotic form of the heat kernel expansion (Q. 

In order to make contact with the standard renormalisation group proce- 
dure that uses an energy scale cut-off note that r c ~ (length) 2 = (energy) ~ 2 
in h = 1 = c units. The proper time cut-off can therefore be replaced with 
an energy cut-off E c using r c = E~ 2 - The divergent part of La becomes 

divp(L A ) = J rf 4 x{iL c 4 trL + L 2 trLx + trL 2 InL 2 }. (12) 

Previous work— used dimensional regularisation 18 with the space- 
time dimension taken as n = 4 + e with e — > understood. In this case the 
lower limit on the proper time integration can be kept as r = and it is 
found that the divergent part of the effective action La contains a simple 
pole as e — )> given by 

divp(L A ) = —V / ^trL 2 . (13) 
167i^e J 

Comparison of the dimensional regularisation result ffl~3l with that of the 
cut-off method ([12"]) shows that the coefficient of the simple pole at e = in di- 
mensional regularisation is the same as the coefficient of InL" 1 . The quartic 
divergence (proportional to L 4 ), and the quadratic divergence (proportional 
to L 2 ) do not appear in dimensional regularisation; they are both regulated 
to zero. If the cut-off energy E c is regarded as a fundamental scale in the ef- 
fective field theory, then neglect of these terms could be significant. For the 
consideration of the gauge coupling constant renormalisation it is possible 
to show by calculating the E coefficient that there can be no contribution 
to the charge renormalisation from the term involving E^ in (I12p. However 
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there is a potential contribution from the quadratic divergence (middle term 
of ( 112]) ): in fact, it will be demonstrated that such a divergence is present, 
as found in the original calculation—, and that it alters the result found by 
using dimensional regularisation. 

It is clear from (|P2|) that the central object of importance for deciding 
whether or not quadratic divergences are present lies in the expression trEi. 
Although tr£i is known for special operators^ 1 ^ 1 ^ it is not known for the 
general operators that arise in the present calculation. The full details are 
technically involved and will be presented elsewhere. If dimensional regular- 
isation is chosen, with n = 4 + e, then G l j(x, x) has a simple poleM as e — >■ 
whose residue involves the heat kernel coefficient E\\ 

4ivp{G l j (x,x)) = --^E 1 i j (x). (14) 

It should be noted that dimensional regularisation is only used as a technical 
device for calculating the heat kernel coefficient Ei, and that this is distinct 
from any choice of regularisation employed for the effective action. The 
calculation of E\ has been checked using another method that does not use 
dimensional regularisation. 

The heart of the calculation now involves finding the pole part of G l j(x, x) 
and then reading off the expression for E\. To accomplish this the local 
momentum space approach 38 is adopted that utilises a normal coordinate 
expansion of the operator A\,-. The general form of A i j is 

= (A^fjdadp + (B a y 3 d a + (Cy, (15) 

for coefficients (v4 a/3 )\,-, {B a ) l j and (C) l j that depend on the spacetime coor- 
dinates through the background field. Normal coordinates are introduced at 
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x' with x^ = x' M + and all of the coefficients in ( II 5p are expanded about 
1/^ = 0. This gives 

oo 

(a^Yj = (AfYj + X^V-tJV 1 ■ • • ( 16 ) 

n=l 

with similar expansions for {B a ) l j and (C)Y The Green function is Fourier 
expanded as usual, 

&iM = J J^n^Mi (17) 

except that the Fourier coefficient G l j(p) can also have a dependence on the 
origin of the coordinate system x' that is not indicated explicitly. 

By substituting ( TlTl) into (J7]) and using the expansion ( fTBT) and similar 
ones for B a and C it is possible to solve for G l j{p) as an asymptotic series 
in 1/p beginning at order p~ 2 . If 

G>) = Go^ip) + GSjip) + GV,(p) + • • • (18) 

where G r l j (p) is of order p~ 2 ~ r as p — > oo it is easy to see that to calculate 
the pole part of G% ( 4 only terms up to and including G^jip) 

are needed. These can be found iteratively beginning with Go l j(p). The 
evaluation of these terms is extremely complicated and most of the brute 
force calculation was done using Cadabra 39 . 

The net result of this lengthy calculation is that the gravity and gauge 
field contributions result in 

2 /3 3 1 2 3, 1 , 3 , 1 ^ 1 2 a - 2 

XxE x = k ( uj + -w 2 + ~£ - + -ujC uEC H u 2 C F 

1 V8 4 8 8^ 2 s 8 S 32 ^ 32 V 

+ (12 + 8£ 2 + 3v + < 2 )A. (19) 
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One check on this result is that all terms that involve l/£ that arise at 
intermediate stages of the calculation cancel from the final result. Another 
useful check is that if the choice £ = 1 = ( is made then the result for E\ 
follows from a standard result— and agrees with the result for trEi found 
above. It is noteworthy that the parameter v that marks the difference 
between the standard and gauge condition independent effective actions does 
not enter the term in the E\ coefficient that multiplies the field strength F 2 , 
and therefore cannot contribute to the quadratic divergences responsible for 
charge renormalisation. This was noted earlier-22 using a different approach. 
(The parameter v does occur in the term involving the cosmological constant, 
but this can play no role in charge renormalisation.) To get the result of the 
standard background-field method, that must also agree with an analysis 
using Feynman diagrams and normal Feynman rules, we simply set v = 0. 
The coefficient of F 2 in this term can be seen to depend explicitly on the 
two gauge parameters £ and (, as well as on the gauge condition parameter 
co. In this gauge condition dependent case it still follows that there is a 
quadratic divergence, a result that is at variance with using a momentum 
space cut-off—. 

A similar procedure can be applied to the ghost field contribution (second 
term of and the Dirac field contribution (last term of ([2}). The Dirac field 
has no quadratic divergence since^ trEi = 0, but does have a logarithmic 
divergence coming from tiE 2 . The ghost field operator results in 

tiEt = ~^ 2 F 2 . (20) 
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Divergences and renormalisation group 



The overall result for the quadratically divergent part of the complete one- 
loop effective action ( TTTT) that involves F 2 is 

_,(!) k 2 K 2 /3 3 5 2 3, 1 , 3 , 1 _ 1 2 a /" , 4 - 2 

(21) 

This conclusively demonstrates the gauge parameter and gauge condition 
dependence of the standard result. If £ — >■ 0, ( — > 0, u — > 1 are taken to 
obtain the gauge condition independent result as discussed above, the non- 
zero result 

r (i) _ k 2E c f >4 p2 (22) 
1 ^ ad " 128^ J d Xt {ZZ) 

is found. 

This is not the complete divergent part of the effective action that involves 
F 2 because there is still the logarithmic contribution. It is possible to extend 
the method of calculation described above to calculate E2, but the results 
can be deduced from the known— >^ poles found for the first two terms of 
(jlip . and for the Dirac contribution 34 . The net result for the divergent part 
of the effective action that involves F 2 and therefore contributes to charge 
renormalisation is 

di ^< r,i, » - 1§ in * + - e <) I < 23 > 

From this the renormalisation group function in ([T]) that governs the running 
electric charge to be calculated to be 

^• e > = T^-3^ (B2 + f A)e - (24 > 
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The first term on the right hand side of (j24j) is that present in the absence of 
gravity (found by letting k — > 0) and results in the electric charge increasing 
with the energy. The second term on the right hand side of (jUj) represents 
the correction due to quantum gravity. For pure gravity with no cosmolog- 
ical constant, or for small A as present observational evidence suggests^, 
the quantum gravity contribution to the renormalisation group /3-function is 
negative and therefore tends to result in asymptotic freedom, in agreement 
with the original calculation—. 

Outlook 

Although the calculation has been done for quantum electrodynamics, simi- 
lar conclusions follow for Yang-Mills fields with and without further matter 
present. It is possible that the running of the gauge coupling constants (found 
by solving ([1]) for a realistic gauge theory) can lower the unification scale in 
comparison to what is expected in the absence of gravity—. In addition, if we 
allow spacetime to have more than four dimensions, then the presence of ex- 
tra dimensions could lower the effective gravitational scale below the Planck 
scaled and render the results of running charges measurable 1 ^ in the LHC. 
There is also an intriguing connection with the weak gravity conjecture^ 
which predicts that the natural cut-off should be a couple orders of magni- 
tude below the Planck scale. The weak gravity conjecture translates into the 
requirement that the gravitational contribution to the renormalisation group 
/3-function should be less than that found in the absence of gravity 43 . 
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